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CKM global fit and constraints on New Physics in the B meson mixing 

O. Deschamps [CKMFitter group] 

Laboratoire de Physique Corpusculaire, IN2P3/CNRS, Clermont-Ferrand, France 



An update profile of the CKM matrix is given, providing numerical and graphical constraints on the CKM parameters 
in the Standard Model. Constraints on additional parameters accounting for possible new physics contribution in a 
model-independent analysis are also reported with emphasis on the leptonic decay of the B d and on the B s mixing 
phase. 
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1. INTRODUCTION 

Within the Standard Model (SM), the quark flavor mixing is described by the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix [1]. The 3x3 unitary CKM matrix can be fully parametrized by four independent parameters among which a 
single non-vanishing complex phase accounts for the violation of the CP symmetry. 

Inspired from the one proposed by Wolfenstein Q the following parametrization, phase-convention independent and 
unitary-exact to all orders in A [3j, is used throughout this document: 
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While the parameters A and A are accurately determined through the measurement of \V u d\ from the super-allowed 
nuclear transitions, \V US \ from the semileptonic kaon decay and \V c b\ from the semileptonic B decay with charm, 
the (p + iff) complex parameter, being the apex coordinates of the unitarity triangle (UT) related to the first and 
third quark families, is less constrained. The metrology of the UT apex via the determination of its sides and angles 
allows to measure the size of the CP violation and to validate the overall consistency of the KM scenario within the 
Standard Model. Any inconsistency would suggest contributions from physics beyond the Standard Model. 



2. THE CKM FIT INPUTS 

s The CKM global fit is performed within a frequentist statistical approach including a specific treatment to handle 
the theoretical uncertainties on some of the inputs (RFit Q). The inputs used for the fit are the most recent available 



input 


source 


value 


(ref) 


update 


UT constraint 


associated theory [parameters] 


\v ud \ 


Nuclear decays 


0.97418 ± 0.00026 


(f4]) 




\v ud \ 


[l^ud form factor] 


\Vus\ 


SL kaon decays 


0.2246 ±0.0012 


([5]) 




\ 11 s 


[|Vus| form factor] 


\Vd,\ 


SL charmed B decays 


(40.60 ± 0.35 ±0.58)10 -3 


([6]) 




\Vcb\ 


[| V^t,| form factor and/or OPE] 


\v ub \ 


SL charmless B decays 


(3.76 ± 0.10 ±0.47)10- 3 


([6]) 




\v ub \ 


l\Vub\ form factor and/or OPE] 


BR(B+ -> t+u) 


Leptonic B decays 


(1.73±0.35)10~ 4 


W) 


* 


\v ub \ 


leptonic amplitude [/s a , /s a / /s d ] 


Am, 


B S B S mixing 


(17.77 ±0.12)ps- 1 


([8]) 




\VuV* b \ 


AB=2 amplitude [B s , /s s , fht , rjg] 


Am d 


B d B d mixing 


(0.507 ±0.005)ps- 1 


(f6]) 






AB=2 amplitude [B B JB Bd f B Jf Bd , Vb] 


\*k\ 


KK mixing 


(2.229 ±0.010) 


(f9]) 


* 




AS=2 amplitude \B K , VB , Vet, Vtt] 


sin(2/3) 


Charmonium B decays 


(0.672 ± 0.024) 


W) 


* 







BR & A CP 


B — * 7T7T, pp, pir decays 


B-factories average 


([10]) 




a 


SU(2) 


BR & A C p 


B — > DK decays 


B-factories average 


([10]) 


* 


7 


GLW/ADS/GGSZ 



Table I: Most relevant inputs for the global CKM fit. The left part addresses the observable inputs with their value and 
reference. The right part addresses the constraint on UT parameters derived from the observable. The upper part of the table 
contains CP-conserving observables while the CP-violating are listed in the lower part. 
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value of the observables from the K and B meson physics for which the theoretical conversion into CKM parameters 
is under control. The most relevant input observables of the global CKM fit, with their value and reference, are listed 
in the table HI The observable values that have been updated for this analysis with respect to the previous similar 
analysis 13| are indicated. 



The right side of the table addresses the 
associated theory, possibly including phe- 
nomenological parameters, needed to derive 
a constraint on the UT parameter(s) from 
the input observable. The value coming 
from the lattice QCD calculation (LQCD) 
for some of the phenomenological quantities 
used in this analysis are listed in the second 
column of the table [ill For comparison pur- 
pose, some alternative and more recent values proposed in the reference [14j are quoted in the last column (note that 
the statistical and theoretical uncertainties have been assumed quadratically summable by the authors). 
A more detailed review of the CKM global fit inputs can be found in [3]. 



observables 


values used for this analysis [H, ,12] 


possible update [14] 




(268 ± 17 ± 20) MeV 


(245 ± 25) MeV 


Ib s 1 !b a 


1.20 ±0.02 ±0.05 


1.21 ±0.04 


Bb s 


1.29 ±0.05 ±0.08 


1.22 ±0.12 


Bb 3 1 B Bd 


1.00 ±0.02 


1.00 ±0.03 


B K 


0.78 ±0.02 ±0.09 


0.75 ±0.07 



Table II: Phenomenological quantities obtained from the LQCD calculation. 



3. THE SM CKM FIT RESULTS 

The left part of the figure [1] displays the result of the global fit together with the 95% CL contours of the individual 
constraints. 



0.0 




Figure 1: Left figure: individual and global constraints in the (p,fj) plane from the global CKM fit. The hashed region of the 
global combination corresponds to 68% CL. The constraints from the experimentally dominated observables and observables 
inducing larger theoretical uncertainties are shown on the right-top and the right-bottom figure, respectively. 

A satisfactory agreement is observed from the various individual contributions at the 2a level, establishing the KM 
mechanism as the dominant source of the CP violation in the B meson system. 

As shown on the right side of the figure [TJ a slight tension is however revealed when comparing the global constraint 
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coming from the observables dominated by the experimental measurements (UT angles) with the constraint derived 
from the observables inducing larger theoretical uncertainties (|V^b|, BR(_B+ — > t + v), \ck\, Amd/Am s ). This 
classification of the observables is for illustration purpose and should not be used to draw any conclusion about the 
possible origin of the tension: a similar discrepancy is also visible when comparing the constraint from the tree-level 
processes with the mixing loop-induced or the constraint from the CP-conserving observables versus the CP-violating 
ones. 




0.5 0.6 0.7 0.8 0.9 1.0 ' 0.0 0.5 1.0 1.5 2.0 2.5 ' D.03 0.04 0.05 0.06 0.07 0.08 0.09 0.10 



sin 2(3 br(b^tv)xio 4 BR ( D ,->w) 

Figure 2: Left: Confidence Level contour in the 2D (sin(2/3),i?7?(i? + — > t + v)) plane for the fit prediction. The cross indicates 
the 68% CL interval for the direct experimental measurement. Middle: the corresponding ID projection for BR(B + — » t + v) 
compared to the direct measurement. For comparison purpose, the prediction for the branching ratio of the D s —> t + v leptonic 
decay is displayed on the right plot (the decay constant value fo s = 241 ± 3 MeV has been used [ljj). 



As illustrated on the left plot of the figure [2J the bulk of the tension is located in the correlations between the CP- 
conserving, UT-side related, theory-dependent observable, BR(i? + — > t+v) (central value of which slightly increases 
with the summer 2008 update 00 an d the CP-violating, UT-angle related, theory-free observable, sin(2/3) (central 
value of which slightly decreases [6j). 

Quantitatively, the minimal \ 2 of the global fit decreases by 2.9 a (2.6 a) when removing BR(i? + — * t+v) (sin(2/3)) 
from the list of the fit inputs. Obviously, the disagreement is not large enough to exclude the possibility of a statistical 
origin of the tension. However, it is worth mentioning that a similar discrepant pattern is also observed for the D s 



leptonic decays (see for instance 20l.l2ll|). suggesting a common origin. As an illustration, the expected and measured 
branching fraction of the leptonic B+ and D s decays to t+v final state are displayed on the right side of the figure [2l 
The possibility of a generic problem with the lattice prediction for the theoretical parameters related to the B meson 
has been investigated. Although the tree-level expression for the amplitude of the the B + — > t+v annihilation decay 
is directly proportional to the decay constant fg d , it can be shown that an under-estimation of the lattice prediction 
for the parameters product fg s x (fB d /fB,) as used in this analysis, can not by itself explain the observed tension. 
Indeed, the BR(_B+ — > t+v) together with the Am^ observables provide a /^-independent constraint in the (p,fj) 
plane in which a clear tension remains as illustrated on the figure [3] for the two set of the LQCD parameters values 
quoted in the table [TTJ 
More precisely, the /s-independent ratio: 

BR(B+ -> t+v) _ 3tt m 2 T ( m 2 T \ 2 /sin/3 



Am d 4 ml v S{xt) V* rn 2 B+ ) Tb+ J \V u d\ 2 B i _ 

allows for a theory parameter-free prediction of the bag parameter, Bb a , from the experimental measurement of 
j3, 7 (or a = tt — j3 — 7), \V u d\, BR(B + — » t + v) and Amj. The resulting prediction is 2.4 a away from the 
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Figure 3: Constraint in the (p,fj) plane from the BR(B + — > t + v) and Airij observables only using the LQCD parameters 
value quoted in the 2nd column of the table [TT] (left) and using the alternative values quoted in the third column (middle). In 
the latter case, both the theoretical and statistical uncertainties are assumed to be Gaussianly distributed, resulting in a more 
aggressive impact of the LQCD parameters. On both plots the global constraint from the whole set of observables is indicated 
by the shaded area around the apex of the unitary triangle. Right: theory-free prediction for the bag parameter Bs d (see 
text) compared to the LQCD calculation quoted in the second column of the table ITT1 



lattice calculation [ill, ll2| as shown in the right plot of the figure [3J This deviation is essentially dominated by the 
experimental uncertainties on the B + — > t + v branching ratio and on the UT angles a and 7. 

A lattice origin of the tension would thus involve the bag parameter Bs d that controls the correlations between 
BR(£? + — > t + v) and the angle (3, as explicitely shown in the above formula. Further investigations are needed to 
check whether the possible correlations between the theoretical errors affecting the lattice parameters, generally not 
provided by the lattician community and therefore neglected in this analysis, could account for the observed tension. 
Eventually, the hypothesis of new contributions beyond the SM is discussed in the next section. 



4. CONSTRAINTS ON NEW PHYSICS IN THE B d , MESONS MIXING 



4.1. Fit procedure 



New Physics (NP) is expected to affect the amplitude of the neutral B mesons mixing in many scenarii (see for 



r-i 



instance [16j]). Assuming NP to contribute mostly to the short-distance part of the AF = 2 processes, a model- 
independent parametrization has been proposed 



17 



18|: 



(B q \M?r NP \B q 



A™(B q \M™\B q ) 



where the label q stands for the d or s flavor of the neutral B meson and the complex parameter A~ p = |A~ p |e « 
accounts for the NP contribution. Assuming in addition that the tree-level mediated decays proceeding through a Four 

Kyi, 7 and 7(a) = TT-Pcs-a 1 



15 



22j), the observables \ Vij 



Flavor Change get only SM contributions (SM4FC hypothesis 
are not affected by the NP contribution and can be used in a (SM+NP) global fit to fix the SM parameters. The 
oscillations parameters, the weak phases, the semileptonic asymmetries and the B meson lifetime-differences are 



1 where /3 C c means the /3 UT angle extracted from the analysis of the charmonium meson decays. 
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affected by the phase and/or the amplitude of the NP contribution as quoted in the table Hill and allow to constrain 
the NP deviation to SM parametrized with A^ p . 



parameter 


prediction in the presence of NP 


Am 


IA NP I v A™ SM 


2(3 




2(3 S 


20?* - 


2a 


2 ( 7r -^M_ 7) _ $ NP 


* 12 „ = Arg[-^] 
Asl 


$?2 m , + < p 


M HM x | A NP 




2|r 12 1 , 2 Jxcos($i + ^ p ) 



Table III: Theoretical prediction for the B physics observables in the presence of NP in mixing. Note the opposite sign 
convention for the Bd and B 3 mixing phase. 



The resulting constraints in the (Re(A^ p ),Im(A^ p )) plane are summarized in the left and right side of the figured] 
for the Bd and B s case, respectively. For a more complete review of the constraints from the global fit of NP in the 
B meson mixing see 231 ] ■ 




Figure 4: Individual and global constraints in the (Re(Ag P ),Im(A, p )) plane for the (SM+NP) fit. Due to the large uncer- 
tainties, only the 68% CL are shown for the individual constraints. Both 68% and 95% contours are shown for the global 
constraint. 



4.2. The B d case 

For the Bd meson case, the dominant constraints come from the mixing phase, /3, and the oscillation frequency Am<j. 
Both reasonably agree with their SM prediction. A 2.1 a deviation is obtained for the 2-dimensional SM hypothesis 
AJJ P = 1. The 1-dimcnsion hypothesis based on the phase only, $^ p =0, corresponding the phase expectation in 
the Standard Model as well as in the Minimal Flavour Violation scenarii, results in a 1.5 a deviation. This slight 
disagreement with the SM hypothesis is essentially the translation of the tension in the SM global fit discussed in the 
previous section (obviously the discrepancy in the D s leptonic decays mentionncd in the previous section is unaffected 
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by the NP scenario in the Bd mixing assumed here). 

When removing BR(£> + — ► t + v) observable from the fit, the deviation reduces down to 0.9 a for both the ID null 
NP-phase hypothesis and the 2D SM hypothesis. 

4.3. The B s case 

The fit for New Physics in the B s mixing also exhibits a 2.1 a deviation for the SM 2D hypothesis, A^ p = 1 and 
a 2.5 a for the ID null NP-phase hypothesis, <&" p =0 2 . Contrarily to the Bd case, the source of the discrepancy is 
obvious: the deviation is fully dominated by the direct measurement by the CDF and the DO experiments of the 
correlated weak mixing phase and lifetime difference, (AT S , 2f3 s ) through the time-dependent angular analysis of the 



B s — > J/>I'<f> decay . Both collaborations report a large phase with respect to the SM expectation [26[. The Heavy 
Flavor Averaging Group [6[ obtains a 2.2 a deviation for the SM hypothesis in the DO/CDF (Ar s , 2j3 s ) combination. 
The other contributions to the NP fit only provide weak constraint: 

• the oscillation frequency, Am s , only sensitive to the module of A" p is consistent with its SM prediction. 

• the semileptonic asymmetries, A S SL , and the Flavor-Specific proper-time, Tp S do not display any significant 
sensitivity to the NP fit parameters in the current state of their determination. 

• the theoretical lifetime difference proportional to the cosine of the weak phase: Ar s ~ 2|Tf 2 1 x cos($" p ) slightly 
tends to push the NP phase toward the zero value as the measured Ar s is larger than the SM expectation for 

2|rf 2 | 



5. SUMMARY AND CONCLUSIONS 

The current data successfully fit the KM mechanism establishing it as the dominant source of CP violation in 
the Bd sector. The current limitation of the CKM consistency tests are of three origins: LQCD quantities, the UT 
angle 7 and the B s measurements. With the recent update of the inputs, a slight tension appears in the global 
CKM fit within the SM hypothesis. This tension is mainly driven by the larger BR(£> + — > t + v) updated value and 
the smaller sin(2/3) value that B-factories recently reported and results from a non-trivial correlation that we have 
explicited in the ratio BR(£> + — > tv)I A TO<4 . Both updates are consistent with the previous measurements but slightly 
go away from their SM prediction. Besides the possible statistical fluctuation of the experimental measurements, 
further investigations are necessary to understand the role LQCD quantities may play in the origin of the tension. 

Assuming the presence of New Physics in the B mixing, a model-independent fit procedure has been applied. As 
a translation of the observed tension in the SM CKM fit, the SM hypothesis is 2.1 a away from the preferred NP 
fit solution in the Bd case. A 2.1 a deviation to SM is also observed in the B s sector, fully dominated by the direct 
experimental measurement of the B s mixing phase by the Tevatron experiments. 

Within a joint fit for NP in the Bd and B s mixing, the 4-dimension SM hypothesis (A^ p = A^ p = 1) results in a 
2.9a deviation with respect to the prefered fit solution. 



A more complete update with all the data eventually available at the end of the summer 2008 is to appear in 15 1. 
In the near future, significant improvements in the lattice QCD predictions as well as more precise experimental 
measurements can be expected in both the Bd and the B s sectors. These data will be scrutinized accurately to 
confirm or not the appearing tension in the Bd sector and the large mixing phase observed in the B s mixing. 



2 A deviation in excess of 3 u in the fit for NP in the B s mixing had previously been reported by the UTFit collaboration ]2-l| and has 
recently been updated to a lower value 12 TJI 
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Eventually a major step in the consistency tests of the KM mechanism will arise with the imminent LHC era and in 
particular with the dedicated LHCb experiment. 
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